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We show that very large negative magnetoresistance can be obtained in magnetic trilayers in 
a current-in-plane geometry owing to the existence of crossed Andreev reflection. This spin-valve 
consists of a thin superconducting film sandwiched between two ferromagnetic layers whose mag- 
netization is allowed to be either parallelly or antiparallelly aligned. For a suitable choice of struc- 
ture parameters and nearly fully spin-polarized ferromagnets the magnetoresistance can exceed 
—80%. Our results are relevant for the design and implementation of spintronic devices exploiting 
ferromagnet-superconductor structures. 

PACS numbers: 74.45.+c,72.25.-b,85.75.-d 



Giant Magneto Resistance (GMR) is the pronounced 
response in the resistance of magnetic multilayers to an 
applied magnetic field 0, S S H ■ This phenomenon 
has prompted a very large interest owing to its broad 
range of applications, spanning from magnetic recording 
to position sensor technology, and to the fundamental 
interest in spin-dependent effects A magnetic multi- 
layer consists of an alternating sequence of ferromagnetic 
(F) and non-magnetic layers (N). The relative orienta- 
tion of magnetic moments in the F layers can be driven 
from antiparallel (AP), in the absence of external field, 
to parallel (P), with a small (up to some hundreds of Oe) 
magnetic field. GMR was originally demonstrated |fjj 
in Fe/Cr multilayers with current flowing parallel to the 
planes, the so-called current-in-plane (CIP) configura- 
tion. In the CIP measurement the magneto-resistance 
(MR) ratio, defined as the maximum relative change in 
resistance resulting from applying the external field, is 
typically around 10% for a number of layers of the or- 
der of 50 - 100 @. These values can be increased up to 
~ 100% in the case of current flow perpendicular to the 
multilayer plane (CPP configuration) 6]. 

In this Letter we show that the limitations of the CIP 
configuration can be overcome by employing a supercon- 
ductor (S) in the non-magnetic portion of the multilayer. 
The use of superconductors in spintronics is not new. As 
a matter of fact, superconductors were used already in 
the very first CPP experiment [|J in order to minimize 
the extra resistance introduced in contacting the multi- 
layered structure to the measuring apparatus. The pe- 
culiar properties of FS structures have been studied for 
several years and this field has been recently reviewed in 
Ref. 0. 

The structure we envision (see Fig. 1) consists of two 
identical diffusive ferromagnetic layers (Fi and F2), of 
thickness t-p, separated by a (s-wave) superconducting 
layer of thickness is- The layers are assumed to be in 
good metallic contact and have length L and width w. 




FIG. 1: (color online) Sketch of the FSF spin-valve. A thin su- 
perconducting film is sandwiched between two identical ferro- 
magnetic layers whose magnetizations (yellow arrows) can be 
aligned both in the parallel (P) and antiparallel (AP) configu- 
ration. An electric current (white dashed arrows) is allowed to 
flow through the system parallel to the layers. The schematic 
representation of the spin-valve effect for half-metallic ferro- 
magnets, showing the diagrams of the superconducting den- 
sity of states, is displayed in (a) and (b). (a) In the P align- 
ment, the lack of quasiparticles with opposite spin hinders the 
condensation of two electrons injected from the ferromagnets 
in a Cooper pair in S. As a consequence, the electric transport 
is confined within the F layers, (b) In the AP configuration, 
two electrons with opposite spin injected from the F layers 
can form a Cooper pair within the superconductor thanks 
to crossed Andreev reflection, thus "shunting" the current 
through the whole structure (see text). 



The magnetization of the two ferromagnets is allowed to 
be aligned cither in a parallel or an antiparallel configura- 
tion |8( . The trilayer is connected to ferromagnetic leads 
separated by an insulating layer (light-yellow regions in 
Fig. 1) of the same thickness as the S layer. The magne- 
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tization of the upper F leads is equal to the one relative 
to layer Fi, and analogously for the lower F leads. In 
the CIP configuration, charge transport in the system 
is dominated by crossed Andreev reflection (CAR) lead- 
ing to a dramatic enhancement of the magnetoresistance. 
CAR was analyzed in several papers notably in re- 
lation to quantum information processing [To| | , an d very 
recently it was observed experimentally in FS |ll| and in 
NS |l2j structures. Here we emphasize its potential for 
spintronics. 



Let us first describe qualitatively the principle of oper- 
ation of the present spin- valve. For the sake of clarity, let 



us first consider a half-metallic (i.e., with only one spin 
species) ferromagnet in good metallic contact with a 
S layer. Quasiparticles with energy below the supercon- 
ductor gap can be transferred into the superconductor as 
Cooper pairs only through an Andreev reflection (AR) 
process [l4j. The latter consists of a coherent scatter- 
ing event in which a spin-up(down) electron-like quasi- 
particle, originating from the F layer, is retroreflected at 
the interface with the superconductor as a spin-down(up) 
hole-like quasiparticle into the ferromagnet. Since only 
quasiparticles (electron- and hole-like) of one spin type 
exist in the ferromagnet, no current can flow between 
the F and S layers . Similarly, in the case of the FSF 
trilayer in the P configuration (see Fig. 1(a)), the two 
F layers cannot transfer charge into the superconductor. 
Current is confined to the F layers and it consists of fully- 
polarized quasiparticles. If the S layer is thin enough 
quasiparticles can also tunnel through it (this will oc- 
cur for t$ values up to some superconductor coherence 
lengths (£o))- In the AP configuration (see Fig. 1(b)), 
each of the two F layers can contribute separately to the 
quasiparticle current through the structure just like in 
the P configuration. More importantly CAR does take 
place. In this case a Cooper pair is formed in the su- 
perconductor by a spin-up electron originating from the 
Fi layer and a spin-down electron from the F2 layer. In 
the AR language, this can be described as the transmis- 
sion of a spin-up electron-like quasiparticle from one of 
the F layers to a spin-down hole-like quasiparticle in the 
other F layer. This is now possible since the quasiparti- 
cles involved belong to the majority spin species in each 
of the two layers. A charge current can therefore flow 
through the S layer as a supercurrent, thereby shunting 
the conduction channels in the ferromagnets [lfiq . This 
contribution to the current will dominate at least when 
the structure is long enough and the quasiparticle con- 
tribution in the F layers becomes negligible (note that 
the conductance of each F layer in the diffusive regime 
is proportional to l/L, where I <C L is the mean free 
path). As a result, one can expect the conductance Gap 
of the AP configuration to be much larger than the con- 
ductance Gp of the P configuration. This can give rise 



to a large, negative value of the MR ratio, defined as: 

Gp — Gap 



MR = 
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A simple expression for the MR ratio for half-metallic 
ferromagnets in the diffusive regime can be derived as 
follows. In the P configuration, the conductance is ap- 
proximately given by [l£j 
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i.e., it is proportional to the number N-\ of open chan- 
nels for spin-up electrons of each F layer, and inversely 
proportional to L. In the AP configuration, the con- 
ductance can be roughly separated in two contributions. 
One (G*), due to CAR, is virtually independent of L. 
The other comes from the direct transmission of quasi- 
particles (proportional to (2e 2 /h)(£/ L)N^): 
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with a being a numerical factor ~ 1. As a result, 

MR~l-a-G*A- — , 

2e 2 t 7V T ' 
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negative and large for L ^> I. This is in contrast to what 
expected in a FNF trilayer, where the MR value is posi- 
tive [5| since the AP configuration yields a reduction of 
the structure conductance. For non half-metallic ferro- 
magnets, the charge current will still be dominated by 
CAR, but the effect will be reduced. 

This qualitative understanding of the effect can be val- 
idated by a numerical calculation of the conductance, 
which was performed within the Landauer-Buttiker scat- 
tering approach. In the presence of superconductivity, 
the zero-temperature and zero-bias conductance can be 
written as G = Gf + G^ , where 
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is the spin-dependent conductance [H. In Eq. T CT 
(T°) is the spin-dependent normal (Andreev) transmis- 
sion probability for quasi-particles injected from the left 
lead and arriving on the right lead, while TL a & is the An- 
dreev reflection prob ability for quasi-particles injected 
from the left lead 13 ■ Similarly, T°' and 7^' arc the An- 
dreev scattering probabilities for quasiparticles injected 
from the right lead, e is the electron charge and h is the 
Planck constant. The scattering amplitudes were eval- 
uated numerically by making use of a recursive Green's 
function technique based on a tight-binding version [2flj ] 
of the Bogoliubov-de Gennes equations 



n 

A* 



A 

-W 



= E 



(6) 



3 




3456789 10 11 0123456789 



yt M£ C 



FIG. 2: (color online) (a) Conductance in the P (black circles) 
and AP (red triangles) configurations versus ts with tp = 
5. (b) Resulting MR ratio. Data were obtained assuming 
L = 150, V = 100%, and U = 8 (see text). In (a) the 
error bars correspond to the standard error over all disorder 
configurations. Lines are guides to the eye. 

where TL is the single-particle Hamiltonian, and u (v) is 
the coherence factor for electron- (hole)-like excitations 
of energy E, measured from the condensate chemical po- 
tential \i. Within the tight-binding description, 7i and 
A are matrices with elements (T~C)ij = tiSij — 7&fi,j} an d 
(A)ij — AiSij, where ti is the on-site energy at site i, 7 is 
the hopping potential and A< is the superconducting gap 
( {...} stand for first-nearest-neighbor sites). In particu- 
lar, ei = £s in the S region, ti = ti in the insulating bar- 
rier, and ti = ep = es T ^cxc in the F layers, h exc denoting 
the ferromagnetic exchange energy, with upper (lower) 
sign referring to majority (minority) spin species. A, is 
assumed to be constant and equal to zero-temperature 
gap (Ao) in the S region, and zero everywhere else. Note 
that this is realistic when the S layer thickness is larger 
than £ H3- Furthermore, disorder due both to impu- 
rities and lattice imperfections is introduced by the An- 
derson model, i.e., by adding to each on-site energy a 
random number chosen in the range [— U/2, U/2], being 
U a fraction of the Fermi energy. In what follows we shall 
indicate energies in units of Ao, and lengths in units of 
the lattice constant a (of the order of the Fermi wave- 
length) . 

In order to analyze the behavior of conductances and 
MR as a function of the various parameters we used a 
two-dimensional (2D) model of the structure, i.e., we as- 
sumed a single lattice site in the z-direction (see Fig. 
1). In our calculations the tight-binding parameters were 
chosen to describe metallic materials: es = 20, t\ = 10 3 , 
7 = 10, so that Co = (2a/7T)^4( 7 /A ) 2 -eo7/A^ = 9.0. 
We set U = 8 and L = 150, so that the F layers are in 
the diffusive regime. To avoid a self-consistent calcula- 
tion of the superconducting gap, we limited our analysis 



FIG. 3: (color online) (a) Conductance in the P (black circles) 
and AP (red triangles) configurations versus tp with is = 40. 
(b) Resulting MR ratio. The same parameters as for Fig. |2] 
were used. 

to values of tg > 30 (corresponding to ~ 3.3£o) dHHH] 
In addition, the conductance was calculated performing 
an ensemble average over 100 realizations of disorder. 

The conductance and MR dependence on S layer thick- 
ness is shown in Fig. [21 Here we chose the ferromagnetic 
thickness tp = 5 and h cxc — 20 (the mean free path turns 
out to be I ~ 21). For this latter value the ferromag- 
net polarization (V) is equal to 100%. Figure |2[a) 
shows that in the P configuration the conductance Gp 
is initially slightly decreasing and roughly constant for 
ts > 5.5£o- This is due to the fact that quasiparticles in 
the two F layers (for large enough ts values > 5.5£o) are 
decoupled, but some direct tunneling can occur through 
thinner S layers. In the AP configuration, the conduc- 
tance Gap decreases until the value ts — 8.5£ is reached, 
and thereafter remains almost constant. Such a behavior 
is expected since, on the one hand, for ts of the order of 
some £0 the conductance is dominated by the supercur- 
rent (mediated by CAR between the Fi and F2 layers). 
On the other hand, by increasing ts , the two F layers tend 
to decouple and the current through the structure is only 
due to quasiparticles flowing separately through them, in- 
dependently of ts- The resulting MR ratio is shown in 
Fig. 12b) and exhibits very large negative values around 
-70% for t s ~ 3.5£ and about -25% for t s ~ 6.5£o- It 
is noteworthy to mention that when the S layer is in the 
normal state (i.e., a FNF trilayer) MR ~ (0.7 ± 1.8)% for 
ts - 4.5£ . 

The role the F layers thickness on the conductance and 
magnetoresistance is analyzed in Fig. |3| for fixed ts = 40 
and V — 100%. Figure OJa) shows that the conductance 
in the P alignment increases linearly with tp according 
to the estimate in Eq. In the AP configuration the 
conductance is again linear in tp with the same slope, 
but it is shifted upwards as compared to Gp. This is in 
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FIG. 4: (color online) MR ratio versus V with tg = 40, and 
£f = 5. The same parameters as for Fig. [5]were used. 



agreement with Eq. the difference Gp — Gap ~ G*. 
As a consequence, the MR ratio (see Fig. EHb)) starts 
from ~ —70% at tp ~ 0.5£o and thereafter decreases by 
increasing the value of ip • 

We finally analyze the behavior of MR, for tg = 40 and 
tp = 5, as a function of the polarization of the F layers. 
Figure 0] shows that the value of the MR ratio remains 
smaller than ~ —30% up to V ~ 87% and then grows to 
larger negative values. Highly spin-polarized ferromag- 
nets are thus required for the effect to be maximized. 
The fluctuations present in the MR('P) curve can be as- 
cribed to opening and closing of conducting channels in 
the F layers as well as to size effects. 

A 3D structure was also considered, allowing the sys- 
tem to extend in the z direction (see Fig. 1). The calcula- 
tions, performed for several values of the structure width 
(w), confirmed qualitatively the overall results found in 
the 2D case. We finally stress the importance of a good 
metallic contact between F and S layers. The presence 
of a barrier at the FS interface would indeed lead to a 
suppression of CAR and therefore of the MR value. 

In conclusion, we have investigated theoretically spin 
transport in a ferromagnet-superconductor-ferromagnet 
trilaycr in the currcnt-in-plane geometry. We showed 
that very large and negative magnetoresistance values 
(exceeding —80%) can be achieved. Such an effect re- 
lies entirely on the existence of crossed Andreev reflec- 
tion. The results presented here are promising in light 
of the implementation of novel-concept magnetoresistive 
devices such as, for instance, spin-switches as well as 
magnetoresistive memory elements. To this end, half- 
metallic ferromagnets such as CrC>2 0, [23, NiMnSb 
H3, Sr 2 FeMo0 6 |3 and La 2 /3Sr 1/3 Mn0 3 |27| appear as 
particularly suitable. Also the Ga/Gai_ K Mn x As mate- 
rial system, exploiting a superconductor in combination 
with heavily-doped ferromagnetic semiconductor layers, 
appears to be a good candidate for the implementation 
of this structure, thanks to Gai _ T Mn, As predicted half- 
metallic nature (for x > 0.125) |28j and to its well- 
developed technology [2flj . 
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